
A 




TRAVELING DENSITY WAVES HT A POSIT TVS COLUMN 

/ 

by 

Richard S. Fridy 
Captain, United States Amy 
and 

Andrew J» Amstrong 
Captain, United States Army 



/ 



Submitted in partial fulfillment of 
the requirements for the degree of 



MSTER OF SCIENCE 
IN 

HIYS ICS 



I 



United States Naval Postgraduate School 
Monterey, California 



19 5 7 



I 




h 




« 



TRAVELING DENSITY V/AVES HI A POSITIVE COLNUN 



* ♦ * * * 



Richard S. Fridy 
and 

Andrevr J. Armstrong 






TRAVELING lENSITY TiTAVES IN A POSITIVE COLUMII 

by 

Richard S. Fridy 
and 

Andrew J. Armstrong 



This work is accepted as fulfilling 
the thesis requirements for the degree of 

MSTER OF SCIENCE 
IN 

PHYSICS 
from the 



United States Naval Postgraduate School 



ABSTRACT 



The field of gaseous oleotronios has become one of burgeoning 
importance in recent years due to the vast strides^ made in the field of 
comraunioations and other related fields# Theories abound, but few 
impregnable mathematical descriptions of the phenomenon of the positive 
column have been advanced# An attempt was made to extend one of' the 
mathematical solutions to a more general case# The resulting equations 
were programmed into a digital computer which produced results in- 
dicating the efficacy of the program, but more data will have to be 
digested in order to establish the reliability of the basic equations# 
The writers wish to express their sincere appreciation to 
Professor Norman L# Oleson for suggesting this problem and for encour- 
aging us and guiding our research. V^e wish to thank Professor H# 
Martinez for his guidance and instruction on the computer portion of 
this thesis— without his help this thesis could not have been brought to 
a successful conclusion# We vrould also like to thank our tvro most 
invaluable assistants—our wives, whose imflagging cheerfulness, faith, 
and understanding made this work possible© 
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CHAPTER I 



INTRODUCTION 

This thesis evolved from the analysis of a paper by Professors M* 
Satosi Watanabe and Norman L. Oleson of the U. S* Naval Postgraduate 
School in vdiioh the phenomenon of traveling density -vraves in a 
positive colxran was investigated. It seems appropriate in the beginning 
to define exactly vihat is meant by some of the terms and expressions that 
vri.ll bo used hereafter. This is done expressly to avoid confusion over 
semantics and not as a gratuitous collective insult to subsequent readers. 

In the simplest instance let us assume that vre have a cylindrical 
glass tube filled vrith a pure rare gas at a pressure of say, 12 mm Hg. 

In terms of every day experience this vrould be considered a relatively 
lo\v pressiirej however, most authorities in the field of gaseous conductors 
refer to this as a high pressure. Sealed in either end of our gas tube 
let us imagine that we have tvro electrodes which conform to the shape of 
the tube* Let each of these electrodes be connected to an external source 
of e.m.f. through an external resistor. Now, it is well established that 
if either the e.m.f. or the resistor is varied, a gloiv discharge can be 
induced to occur in the tube. Generally speaking, the glovr discharge is 
a steady state gaseous discharge vdiich appears at reduced pressures and is 
characterized by relatively low current densities and relatively high 
potential gradients; perhaps it is bettor to say in this latter case, 
relatively high values of electric field to pressure ratios, about 
which more will bo said later. This gaseous discharge has a typical 
st3ruct\iro -i^diich is marked by a rather steep potential gradient at the 
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cathode and operates basically by electron liberation resulting from 
positive ion bombardment at the cathode# Quite simply, this amounts to 
nothing more than having the gas act as a conductor. However, the pheno- 
mena occurring within the gas v/hile it is acting as a conductor cannot be 
so simply explained. In fact so many phenomena occur that it is not 
feasible to discuss all of them within the scope of this thesis. As 
indicated in the opening statement, we shall confine our attention to a 
single phenomenon in the positive colxmin of the glow discharge. In order 
to facilitate our development we shall make some simplifying assumptions 
about the composition of our positive column. We shall assume that it is 
composed of three fluids: 

a. A negative ion "gas" which we shall further assume to be 

composed entirely of electrons 

b. A positive ion "gas" which we shall further assume to be 

composed either of atoms or molecules of the gas vrith 
vriiich the tube is filled, but each of v/hich carries a 
single positive charge 

c. The neutral particles of the gas itself—these particles 

are by far in the majority. 

It is customary in discussing gaseous conductors to refer to ionized 
regions containing high and apprpximately equal concentrations of positive 

(2 ) 

and negative ions as the plasma ^ ^ With the proper environment there can 
be caused to appear in the positive column alternate dark and light bands. 
These bands may move or they may be stationary. In either case they are 
called striations . If they appear to move, they are called moving 



2 



striations; if they do not move, they are called standing striations» 



As a parenthetical matter of interest, it seems vrell to note that moving 
striations are not detectablo by the unaided eye# 

The traveling density 7/aves referred to earlier may or may not be 
moving striations# Certainly ono could be excused if in the first glow 
of understanding, he impetuously ascribed the title of moving striations 
to these traveling density waves; however, sober and mature consideration 
would make such a conclusion doubtful# The important point to remember 
hero is that moving striations are an observable phenomena and not an 
ethereal speculation# 

A mathematical approach by Y/atanabe and Oleson has been undertaken# 

In their approach they attempt to provide a theoretical basis for ex- 
plaining moving striations# However, they assert that certain of their 
approximations may bo in error if applied to actual striations ^^ ^# It is 
the intent hereafter to justify or refute these assumptions# In order to 
do this it is proposed to extend the linear treatment of Vfatanabe and 
Oleson to the non-linear case and thereby arrive at a more general solution* 

► 

A search of the literature v/as instituted to garner the latest, most exact 
data available# Y/here these data ^^re found to be either doubtful or non- 
existent, data wore taken from current experiments being performed at the 
U# S# Ilaval Postgraduate Schoolo The ultimate aim was to substitute these 
data in the pertinent mathematical equations and solve them by moans of a 
digital computer# The solutions then could be compared with the results 
obtained experimentally and a reasonable conclusion could be reached as to 
the validity of the original premises# 
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Since we have assumed that v/e have a fluid filled tube, let us estab- 
lish a coordinate system to describe the particles in our fluid# In this 
system the neutral particles are always fixed in space and only the positive 
and negative particles move vfith respect to our coordinate system# 

If P is the density of any physical quantiby, Q the rate per lonit 
volume at which the quantity is being produced, and 5V the volume occupied 
by the quantity; then. 



Let us define dV' as a small volume of fluid and agree that 5V always 



the same fluid particles# In general the volimie, 4V, will then change with 
time# If we assume that SV is in the form of a rectangular box, then let 
us say that its dimensions axe Sx, 4y, and ^z, such that; 



Now, the x-oomponent of fluid velocity, v^^ nay be different for the op- 
posing faces of our box moving in the x-direotion<> If so, ^x will change 
with time at a rate equal to the difference between these two velocities; 
thusly. 



Qiv 




designates a volume element that moves v/ith the fluid suoh that it contains 







and by a similar development we can say. 



we have ; 



Therefore, 



JiV 

d t 



- t 



dt 

- / ^ 






^3r 



4 £V . 

dt - 






d/F 



Inserting this in the equation above for 

svj^*ev-^SV=Q(V 



dividing out /iV, 
■ 



dt 






bi4jby definition, 

v(7f 



therefore, 

ii h V - vpf e7V= Q 

Yfhich by a relation in vector analysis reduces to: 

This of course is nothing more than the equation of continuity. A more 
rigorous derivation can usually be found in most texts on mechanics# 



See for instance K# R# Symon, Mechanic s (Addison-^^esley Publishing 
Compsiny, Inc#, Cambridge, Massachusetts, 1953), p# 265 ff# 
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For the sake of simplicity we shall assume that all of the negative 
ions are electrons and that the positive ions are formed by the collision 
of electrons with neutral particles and by no other process-- a reasonable 
assumption if we know that our gas is not an electro-negative gas^ \ See 
fig# 1 for a simple graphical representation of this ionization phenomenon# 







If we place a charged particle in an electric field, the particle will 
move under the influence of the Coulomb law of attraction# It seems 
reasonable that if there are neutral particles in the field which are un- 
affected by Coulomb forces, the charged particles moving in the field will 
collide with the neutral ones# This obviously occurs# However, since the 
charged particle is still under the influence of the field, the particle 
7/ill again tend to move after collision toward the electrode opposite in 
sign to the particle# The overall effect is to give the charged particle 
a notion generally in one direction (see fig# 2)# 






- - - ^ - 



-0 



o 



fe O (y \fOf/ 

Since the Coulomb forces operate to move the particle, it would 
appear that the same lex{ acting between oppositely charged particles would 
cause them to recombine# This does happen under proper conditions; i#e#, 
high pressures (of the order of one atmosphere)# Ilov/'ever, at lorv/- pressures 
the collision frequency and collision losses are' low and the relative 
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velocity of the charged particles is higher# Because the velocities are 
higher, collision without recombination takes place since the attracting 
field is not powerful enough to deflect the ions and make them move 
closely together over a large enough distance 
Let us nov/- define: 

n” = negative ion density 
xif = positive ion density 
V = negative ion drift velocity 
v^ = positive ion drift velocity 

z = number of ion pairs created per unit time per one ^ 
electron (constant if the temperature variations 
are very small 

then from the equation of continuity: 



<)t 

<it 



Further consideration of a charged particle in an electric field 
will inevitably lead to the conclusion that the velocity of the ion is a 
function of the field* Experimental evidence shows that the velocity is 
proportional to the field to a first approximation, viz: 

£ 

(see figs* 3 and 4) 
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Figiure 3, Drift rolocity of electrons in argon as a function of B/p, 
R. A. Nielsen, Phys . Rev. 50, 950 (1936) 





Figure 4, The drift velocity of Argon atomic ions in Argon as a function 
of E/p • J. A. Hornbook, Phys . Rev . 84, 615 (1951) 
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The use of vector notation is justified ■when consideration is given to the 
fact that the motion of the ion is directed. We may therefore say: 

vrhere K is a proportionality constant called the mobility— a quantity ■which 
will be more fully discussed later. Considering again figs. 3 and 4, it 
can be seen that the magnitude of the velocity of the electron is much 
greater than that of the positive ion; indeed, in certain insteinces the 
velocity of the electron can be shown to be a thousand times as great as 
that of the positive ion. This fascinating facet of the properties will 
be given more light later in the discussion of arabipolar diffusion. 

Envisioning a cloud of ions (or electrons) in the tube, one is led to 
the inevitable conclusion that the charged particles may exist in varying 
densities. If this be the case, then obviously the ions v/ill flow from 
regions of higher density to the regions of lower densrby by diffusion. 

This will of course contribute a component to the velocity at which the ions 
move, and we can say; 



V, 



OtfPOJfOA/ 



~vn 



where is the concentration gradient with the negative sign indicating 
that the motion occurs in the direction of decreasing concentration. 
Therefore, 



'Diffusion/ 



z -1 



vr\ 



iwhere D is a proportionality constant called the diffusion coefficient for 
the gas. A more adequate treatment of the derivation of this last 
equation can be found in some detail in most texts dealing v^ith the subject 
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of diffusion#* Combining these contributions to Telocity vfe have, 
eqn (2) iT=Kf-iw 

The electric field under -vvhich the particles act is composed of t\*ro 
parts, the impressed electric field and the electric field generated by 
the space charges; such that, 

E'-E,fE' 

where Eq is the impressed electric field and will be considered constant, 
and E* is the electric field generated by the space charges and depends on 
the number of electrons and ions present; i.e#, the density# Therefore 
can say (from Poisson’s Equation): 

- V-E'^ 

v/here e is the electronic charge# 

It is proposed to solve equations (l) by assuming that the actual 
solution is a superposition of a steady uniform solution n^ plus a non- 
steady, nonunifom disturbance 2 /. In order to linearize the equation 
Vfatanable and Oleson assumed that the nonstoady nonuniform disturbance was 
small compared to the uniform, steady disturbance^^ which is not in agree- 
ment vfith experimental results# In the following development no 

such assumption will be made— to the contrary, we shall assume that at 
some particular tinp the orders of magnitude of the steady, uniform and 

See for instance L, B. Loeb, Basic Processes of Gaseous Blectronics 
(University of California Press, Berkeley and Los Angelos, 1955), p.l88 ff. 
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nonsteady, nonuniform parts may be approximately equal# This is not to 
imply that they are always equal, but rather that they could be equal at 
some particular time# The vford ‘^uniform** as used here means that the 
qusintities do not depend on the x-ooordinato taken along the axis of the 
tube^^)# In keeping v/ith our proposed solutions let us nowwite: 



oqn (3) 






n - n^tP 

v/e define £ as the field generated by the steady, uniform charge 
densities n^# 

Substituting the value given for V in equation (2) into equation (l) 

irn^o 

^n'})] ' ^rf-0 



clearing the parentheses: 

+ V P rfJ^F - (v }] - Hn ^ 0 

at ^ 

Substituting equation (3) in the preceding equations; 

^ ^ KB - if KB - P'7/?;- d~Vv ]- 2n;~ iv ro 

7. [r\:\CB t vk"b -D^vn!-DWv^J-zn;' 2 zT -o 

Directing our attention to ‘the terms containing n^ and recalling that by 
definition these are constant; therefore, 

^-0 and 
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vfhich leaves, 

— -f- V (:n:Ke -7JKe-D'v2Jj-2n:-i2j'^o 

— ^ 7- ■*■ ~ ^ 7V^2 - 2 n;- 2 jS ~0 

it ^ 

but v/e have stipulated that E a Eq •/ E’; it then follows : 

^ t7-[^n'K£o~ nll{ e'~ zf KEo~ V l< - P' 7zi J - 2 n:~^ZJ'-o 

V fn^K^E. + n.K X' t t/K % - P^K "X- i>V -^r>:-2p^0 

Again vre refer to the fact that Eq is a consteint, consoquently, 

V-Eo-O 

Omitting the zero terns, our equations become: 

^ ^7'[-nlKE'-2^X'E(,-7f X'e'- D'vv]' ini - 2 pzo 

3 u 

7- f?: kX ' t kXo ^ pX X'-p" 9p"J - z^o'' 2 ^0 

We have defined C as the field generated by the uniform, steady charge 
density n^# Let us now define ^ as the field due to the nonuniform, 
nonsteady disturbance ; such that. 



Inserting this relation into the equations above: 



'^ 7 / 

it 

it 



+ KE,-p~Ki’ vKe-P' 7 2 v-o 

t [n: K 1 1 n! Ke ^ uXX, ^ - ^^=0 
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Removing brackets, 

If *'i;K 7‘i -n;KV- K-pS7-e -Kt-Vv-KT/v-Q-Ke’^v 

- 0' v^iJ- Bn; 

^ t e f r + k\ V' Vt Kv* 7- ttKl^ 7z/rf<W 7.e^K'e-V 2J' 

' =0 

From Poisson’s Equation: 

V-^ ~^Te(n;-ii;) a/?cl V'T Z ‘irre(i)*-7S) 

v/here e is the charge on the electron# Recalling that our region of 
interest is near the axis of the tube removed from the ends and assuming 
cylindrical symmetry, v/e find: 

, 1 , 



V7J^ 






Vfith these conditions in mind let us rewite the equations: 

^ - n; f< [4 TTt {nl-nl )\-nlK re {if- zJ‘)J-k:E^^ - KiSlHrre n; 

-Kif [i ire iv^- w - 2 2^’ =0 

^ t »5,'k'[^ treint-n;)]* n:K'[tre(v’-p^->!^S, ye (.it- O] ^ 

[t ire (u'-- v-f! +K'r. D’ - ?«,--? ^ o 

Inasmuch as our interest in the field is in the central region along the 
axis of the tube, we note that the field t (or more particularly ) will 

bo zero because of the steady, unifomity of the charge density nQ causing 
it# liirther, in this region n^ is approximately equal to n^# 

^ ' 5 ?' ' in: --0 



eqn (4) ^ 



tK's^ ^t\{v-v)v-in:--C> 
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vrhero. 



Y-^Tre. K\zn:-n:)-2 
HrreK n; 

')yHireK\zri:-r^l) 



)(^z Hire t i 
- Hve^ 

\ ' HrreK'' 



For canparison purposes the terms. 




\Thich appear in equations (4) do not appear in the original paper by 
»Vatanabe and Oleson# These, then, are the terms dropped by V/atanabe and 
Oleson as a result of the mathematical expediency employed by them to 
linearize their equations, and as a result of disregarding certain higher 
order terms* 

Equations (4) are the ones for which constants v/'ill be computed for 
a particular set of conditions; e*g*. Argon gas, 12 mm Hg pressure, etc# 
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CHAPTER II 



EVALUATION OF CONSTANTS 

Reiterating our original contention, vre propose to apply equations (4) 
to an experimental situation* In order to do this certain constant terms 
must be evaluated. Other constants are well knovm but vrill be repeated 
here in order that our results may be identically reproduced. In this 
latter category of constants we find; 
e = 4.8 X 10"^° e.s.u. 
k = 1.38 X 10"^^ ergs/ °K 
= 3.1416 

In the former class of constants there will bo a need to loiow the values of 
K’. D", Eq, z, n~, and n^^. Rather than list them arbitrarily, it 

seems more appropriate to derive them individually in order to indicate their 
reliability# 

The particular experimental situation to v/hioh vfo shall apply equations 
(4) is as follovTs: 

A tube 0#5 cm in radius containing pure Argon at a 
pressure of IE mm Hg with currents of from EQma to 
lOOma in our gaseous conductor^ o 

\l±th these conditions the value of the applied electric field can be found# 

From fig# 5 it can be seen that for a pressure-radius product of 6 we have 
13 

in our case an - ratio of approximately 0#5# At a pressure of IE mm Hg 
this vrill yield an of 6 volts/cm# The value of which we obtained for 
our case yields a value of E X 10 statvolts when converted to e#s#u# 

For a description of the method of obtaining these cxirves refer to the 
source cited in connection vfith figure 5. 
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From probe neasuroment s in the plasma within the tube it has been 
experimentally determined that 

i n^ a 7 X 10^^ ions/om^ 




RP 

^rA Hg") 

Figure 5 

( r o n i ^ e d 6a5€6 ) 

It appears rather obvious that a relation exists between mobility^ 
diffusion, and the number of ion pairs created per unit time per one 
electron. Since tho latter two terms have been defined, it remains only 
to define tho first term — mobility. 1.7e have rather arbitrarily stated 
that mobility is the phoport icnality constant which relates the drift 
velocity to the field. However, it seems more important to have a 

**feeling” for the phenomenon of mobility than to have a pedantic Icnovrlodge 

« 

of the word itself. Consider a^ain a charged particle moving under tho 
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influence of a field and within a region occupied by the neutral particles 

of a^(fig. 2). ^^e can say, then, that mobility is the drift velocity of 

}|( 

an ion moving in a unit field. 

In order to utilize the experimental data available it is neces- 
sary to determine values of mobility for the positive ions and the 
electrons. In dealing with positive ion mobility it was found that there 
is a difference in mobilities depending on whether the ions present are 
"atomic” ions or "molecular" ions. It has been shovm that when pressures 
of the order of 6 ram Hg are encountered in Argon, and ions are found 



In searching for a best value for mobility the phrase "reduced 
mobility" was often encountered. Reduced mobility is best described 
as(7)s 




where p is the gas pressure and T is the gas temperatiire in °K, Many 
results are given in terms of reduced pressure, p^, (see for example 
fig, 4). This reduced pressure is defined as 1^1 i 

273p 

Po = 

Veiy often in studying graphically presented results one observes 
that experimenters and theoreticians will plot ion mobility versus 
field to pressure ratio (i.e., —-where the pressxire may be reduced 
pressure or that of the gasj P e,g», figs, 3 and 4, It can be 
shown(2) that p is inversely proportional to the mean free path, 
thusly; | 

P™” X 

If this be so, then: 




\vhere EA is an expression for the energy acquired in a m.f,p. 

This gives the researcher a more vivid and accurate picture of vdxat 
is happening. 
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to be present in the ratio of 0,5 



(3) 



The mobilities of these two ions 



wore found to be 2,6 eind 1,6 respectively^^ Ascimiing that the 

two types of ions are distributed in the ratio indicated throughout the 
positive ool\jmn, a weighted average mobility was foxmd to be: 



K*^“ 4,04 X 10^ om^/stv-seo 

If the tv/o types of ions are assumed to be present in equal concentrations, 
then; 



= 4,38 X 10'^ om^/stv-seo 

As can be seen, the orders of magnitude are the same or at least near 
enough that either term could be used v/ith some degree of accuracy. 

Prom fig. 3 the drift velocity of an electron vms determined to be 
0,36 X 10^ om/seo. To a first approximation. 



K" = X = z 1.8 X lo"^ omVstv- 

E 2 X 10-2 



sec 



The diffusion coefficients, and D~, are most easily found by 
employing the mobility. The relationship between mobility and diffusion 
is derived in nxmierous texts on the subject^^^. The relation is: 



K 

D 



a 

kT 



Where o is the elootronio charge, k is Boltzmann's constant, and T is the 
ion temper attire in °K, The temperature of the ions is approximately equal 
to the temperature of the gas^^*^^. From this equation, then, the 
following values wore determined: 
vr/s 

3,48 om^/seo 
D “ 6,46 X 10^ cm^/sec 

In order to compute the value of z it is necessary to employ ambi- 
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polar diffusion in the positive column results from the previously shown 
fact that the diffusion coefficient for electrons is greater than for ions. 
Assume that as an original condition there exists near the axis of the 
positive column an almost equal number of electrons and positive ions* 

After a short interval the electrons will tend to diffuse tcwrard the walls 
of the tube at a rate greater than the positive ions* This will cause a 
negative charge to be set up on the walls of the tube. The preponderance 
of positive ions remaining near the axis of the tube (duo to their slower 
diffusion) causes an electric field to be established between the axis of 
the tube and the wall* This field in turn will cause a slovring down of the 
electrons moving tovmrd the wall and a speeding up of the ions* It happens 
that wall recombination is preponderant at the pressure we have specified* 
This is borne out by the absence of recombination spectra indicating that 
volume recombination is not important here. In an equilibrium state the 
losses due to recombination are equal to the rates of production* The 
electric field across the positive column is such that electrons arriving 
from the faraday dark space aro accelerated to produce enough electrons 
and ions per unit length of the tube to compensate for losses, thus 
establishing equilibrium' Specifically, in our case we have restricted 
our consideration to the region near the axis of the tube so that radial 
variations neod not be considered* The point to be remembered here is 
that under the action of the radial electric field, an average velocity 
of diffusion is soon reached such that. 



which leads to the fact that tho total ambipolar current is zero* As vre 
have alreac^seen. 
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Also, wo have previously implied (page 13 ) that near the axis of the 
tube ■ 0 and z 0. Using the expressions for velocity and eliminating 
the field be Ween them, it can be shown that^^^; 

_ d4;“ -/ 

^ ~ K" 

« 

v/here is defined as the ambipolar diffusion oooffioient* Further, 

/ 

since K“ is much greater than K ' and is much greater than (vdiere 
and Tp are the electron and ion temperatures respectively), it can be 
shown that ; 



B D” 

Tp K- ® 

Employing kno\ra data (Tq = 12,550 and Tp s 300 ®K) in the approxi- 
mate formula for yields : 

Dq r 145 cmVssc 

From the ambipolar solution to our problem it can be shown that: 



z 



(M2S1Ld 

(R)" “ 



v^here R is the tube radius® Inserting the value for computed above 
and the radius of our tube-- 0#5 cm^ ^ for the conditions vre have chosen, 
we find: 



z - 3#35 X 10^ ion pairs/electron-sec 
Inserting the constants just computed into equations (4) yields: 

^ + 7 . 6 / /.? ~ J025(7f-v)2J ~Z.3‘i5 



^ % 






3. HU p(t ?08 ). 7/ X + y OV X //)" 5 X. 0 X /O'^ (v- tS) V*- 2.3V5 X ' 0 






<JX 
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Little accuracy is sacrificed by disre^jarding those valuos of constants 
which are dominated by the exceedingly largo values of other constants 
which are additive® 

In order to solve this rather formidable differential equation, 
little error is involved if the partial derivatives aro approximated by 
total derivatives and minutely small increments are talcen to be used in 
a digital computer. This may bo more clearly seen later v/here the above 
equations are set up for programming into the computer# Once the equations 
aro written in a ^’useful form” for the computer, we find we must have some 
starting point; i.e#, some initial values for x, etc# at some 

time, t# Since the approach we have taken appears to be a rather unique 
one, there wore no readily available precedents by which we could be 
guided# Harking back to the original intent of this thesis, the answer 
suggested itself. Since \7atanabe and Oleson^^^ assumed a solution to 
their linear equation, it seemed reasonable to take their evaluation and 
apply it to our problem. If this satisfied our equations, then credence 
is lent to the original hypothesis thereby extending the theoretical con- 

I 

tention to actual conditions* Watanabe and Oleson^^^ assunod linear 
solutions : 

■where a7^ and a” are complex amplitudes, k is the ■wave number, and cj is the 
angular frequency. Thus, another problem arose — the feeding of a complex 
number into a digital computer. It is not impossible to oxocute such a 
program, but in view of the physical construction of the computer, it v/ould 
be necessary to sacrifice some of the accuracy of the original program in 
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order to do so* Hon,vevor, since our solution consists of a real part and 
an imaEinary part; that is to say, a real solution and an imaginary 
solution, either is a solution to our problem* So, \re will program only 
the real solution. This reduces our solution to: 



Again wo have introduced terms which must be evaluated* Directing our 
attention to k, the wave number, and referring to the vrork of Watanabe 
and Oleson^^^, we observe that we musb assume their linear coimterparts 
of our equations (4)* This results in no inconsistencies when the 
restrictions placed on our equations are considered* Y\Tatanabe and 01eson*s 
equations are: 




Letting tho time arbitrarily be zero initially these reduce to: 



a a“cos kXQ 
m aTcoS kXo 





where Yz , cind are the same as tho ones previously defined* 
Substituting 




into Y/atanabe and Oleson’s equations and performing tho indicated 
operations will give: 



7 ^ D“ic^ -f )a“ - aV = 0 
■/ iK/Sok / D^k^ -f )a'/ = 0 
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In order for Wo homogeneous equations to have a solution the determinant 
of ooeffioients must be zero* 



-iw-iK'Egk -f D‘k^ ^ \ 

-ioJ-/ iK^S^k -/d4:V \ 

Upon solving this determinant we have, 

[-iu- iK'S^k -f D"k^ -f |^-iw -/ixT^Eok 4 = o 

Recalling that if a complex equation can be written: 

R / il = 0 

where R is the real part and il the imaginary part, then R - 0 and il = 0 
for the equation to be satisfied© 

Hone particularly for our case, 
for R a 0, 

-o/ -/wK7^Qk-uK"Eok7^K“K7^E^^7^ (D"kV ) - Yi ° 

for 1=0 (since i / O), 

-w[n/kVYjJ -K"EQk[D/k2 / Y3] -cd[D“k2 -/ K^^Eok D”k2 ■/ Y J = 0 

Solving this last equation foroO: 

kEo [}c^(k/d“ - - (K" Yj - Y, )] 

[(D-^ ^ D-)k^/ {i, 4- )] 

Eliminating cO beWeen the Wo equations yields, 

eqn (6) (dV 4 I ){dA^ ^ Yj ) - Yjy 

^ k^E§(K7^ 4 K")^(l/lc^ 4 U )(D~k^ 4 Y. ) = q 
pD^k^ 4 S) 4 (D‘k^ ^>3^ 

Introducing new variables, 

eqn (7) ^ = D"k^ 7 ^ Y, and ^ = d4^ 4 
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.2 . 



Solving for k in the first of equations (7), 



k2 = 



^ Y. 



2 

Substituting this value for k in the second equation (7), 

/ D- ^ ^ ^ 

Talcing cosnizance of the fact that D“ ^ we can say: 

1 - '<1 

Referring to the original expression for ^ , 

7 = d 42 / l(j 

and considering the preceding approximation, vre observe: 

d4:^ = 0 

Since D7^ is not zero (see page 18 ), this approximation can be true if sind 
only if, 

k^ i 0 

Inserting the approximate values of and from equation (7), et seq, 
into equation (s) and observing that k^ £ 0, 

If we map equation (6) in the plane, we have approximated very closely 

a point on the extremely complicated curve represented by equation (6); 

2 (l) 

i#Oo, the point v/here k z 0. It can be shown from an elementary 
geometrical relationship that knowing a point and the slope of a curve at 
the point, the equation of the tangent can be vn:'itten# In our case. 
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where, g= ^ and h = KAr 

Yj D 

Substituting the values of tj and ^ as shown in equations (7) for the 
appropriate terns in equation (8), and noting that: 
z 4:7teli^n 

2 

we solve for k : , 

1,2 r 

4‘ YTVit -/lilEoikT 

47Tne 

Recalling the proposed linear solution, v/e direct our attention to the 
complex amplitudes and v/rito: 



a'/ r Ae (with A and ^ real) 



which can be written: 

r A(oos ^ - isin 



a" 



but as has been shovm, 

(-i 60 - iK'E^k 4 4 )a” “ 

From those last two equations we can say. 



z Acos ^ - iAsin ^ s 

a" 

Equating reals and imaginaries. 



D~k^ 4 r, - i(^4 K'hok ) 



Aoos 

Solving for A, 



_ D^k^ 4 and Asin ^ 3 t*? .. , 7^ . . .K-hfi l. - 
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From which. 



ton ^ s 



(J -f K~Sok 

D-k2 / 



Insorting, 



Q 



z Ck^CK^D" - K~p/) - (K~Yj 

[(D' / IT )k^ •/ { 



and 



Simplifying yields, 

/X ^ (K^ 7^ K")kEo 

eqn (9) tan fi - — 

[(D7^ -/ D-)1c 2 y ( jf, V 1i )] 

It can be shcwn^^^ that tan ^ has one maximum bet^veen 0 and ^2. 

Talcing the derivative of tan ^ with respect to k and equating to zero we 
find at this point: 



eqn (lO) k"^ = 






[2Eo(K‘/ / K")J (O'/ /• D’)^ 

Squaring equation (9) and inserting equation (lO) gives: 



(tanj^) 



2 

max 



( k ’^ -/ r')ho 
[Md-/ -f D-)( I 4 I )] 



Calculating this value of 0 from the known values of the constants 
in the preceding equation, we used this as the maximum value of a phase 
factor with which to enter the computer along vrith our previously described 
boundary conditions. The value of ^ obtained as described was found to 
be 1°54*. 
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ClIAPTiiR III 

ASSU1.IPTI0NS AND CRniQUB OF ASSmiPTIOlTS— 

SOME DIFFICULTIES ENCOUNTERED 

The experimental data yield a Langmuir plot vAiich is a straight line 
over t\To oyoles leading to an assumed Maxwellian distribution for electrons^®} 
The assumption that the ions have a Maxvrellian distribution stems from the 
fact that at several millimeters pressure, the temperature of the ions is 
approximately equal to that of the gas molecules The assumption that 
we havo radial ambipolar diffusion vuidor our stated conditions vrithin the 
tube is based on tho results of Karge, Hooks, and Oleson^^^. In this work 
they state that the instantaneous field intensity is due to a very small 
surplus of positive ions or electrons. Their reasoning for this is very 
amply stated in their work. • Since we can say that n” i , it follo\TS 
that vfo can say tho ambipolar curront at the vmlls is approximately zoro. 
Further, 



dn~ 

dt 



dn 

dt 



■/ 



and v” - v^ 
r r 



vjhero we define v“ and as the negative (i.o., electron) ion and 
positive ion radial drift velocities respectively. 

In our preceding work and tho subsequent computer ^vork it was assumed 



that z vras a constant. If z is also evaluated from tho relation 



(4). 






eqn (11 ) z £ 9 X 10' ape" in>\ e7 

using the temperatures observed by Karge, Hooks, and Oleson'' \ it is 
found that z varies vridely \rith temperature as shown bolov:: 



Te 

16,440 

14,660 

8,200 



1.89 X 10^ 
5.3 X 10^ 
2.22 
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Farris notes a variation in z based on temperatures (the data 
being taken from Pupp’s vrork), Y/hen this information is applied by the 
uso of our equation (11), the follOT/ing results are obtained: 

Te 

28,000 • 2.7 X lo"^ 

11,520 1,63 X 10^ 

It should be noted that the validity of the equation for z depends on 

having HaxiTellian distribution of elootrons iihich we have assumed, 

Farris has recently published a short paper on a theory of moving 

striations based on the diffusion theory. His basic assiimption is quoted; 

If it is novf assumed that the small difference beWoen the 
ion and electron concentrations has the same effect axially 
as radially, and that moan values of mobilities and diffusion 
coefficionts may be used, the theory of a stationary positive 
column may be readily generalized to include changes of 
concentration axially and vrith time. Eliminating the ovirrent 
and field from the usual conservation and current density 
equation gives; S - dn/dt - vdiero S is the rate of 

production of ions per unit volume per second and the 
Schottky apbipolar diffusion coefficient given by: 

Dg - kTfKTye. 

The major assumption different from our assumption is that the axial 
drift velocities of the electrons and positive ions are assmed equal 
within the positive column. Using our equations (l) and oliminating the 
field wo obtain: nv = Dg^n, Substituting this in our equation (l) we get 
Farris’ equation (l). However, we question the validity of the assumption 
that the axial drift velocities of all ions in the positive column are 
equal. The interesting point is the variaH;ion of z with electron 
temperature o 

iYora these considerations ive observe that there is a fallacy in 
our assvnaption that z is a constant. True, we imposed the condition that 
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this was valid only for small tomporaturo variations. But at tho time 
we imaginod this variation could easily be on tho ordor of 100 °I1, How, 
it appears that this is not so. Certainly this is a field for further 
inquiry. How this deviation from our original assumption affects o\ir 
problem will bo shown in the next chapter. We have taken znQ to be a 
constant of the value of 2.345 X 10^^. 

It seems noteworthy to point out that thoro is a phase difference 
botvrecn tho maximum degree of ionization and the maximum degree of 
excitation ifdaich has beon observed in some rare gasos, however this was 
not observed by Karge, Hooks, and Oleson^®^. As a point in fact, they 
observed that a maximum in light intensity vras accompaniod by a maximum 
Tq, and a minimum in light intensity was accompanied by a minimum Tg. 

This vrould suggest that there is no phase difference in the case of Argon 
ions in Argon gas at 12 ram Hg prossuro. If we can assumo that equation (11 ) 
holds and further assume that tho maximum light intensity is also evidence 
for maximum excitation, then v/e can say indeed there is no phase differ- 
ence between the phenomena. 

The temperature dependence of mobilities and diffusion coefficients 
was ignored in this thesis. Since the evaluated constants in equation (4) 
vrere taken to be average values it is suspected that temperature 
dependence is of little significance. 

The relationship between mobility and diffusion used in this thesis, 

K - 0 , (pagelS) may be a tenuous one. Von Engel^^^ states that the 

D - kT 

♦ 

Private oaramunication vd.th N, L. Oleson, v/ho vrith A, ’W, Cooper, has 
found that maximijm light intensity and maximum electron temperature are 
not coincident in neon at several mm prossuro. 
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the equation is not applicable to electrons for at least tvro reasons; 

electrons drift with velocities wliich are not proportional to the applied 

field and their mean free path depends on thoir speed. In defense of our 

use of this equation vfo relied on experimontal data (see fig. 3) which 

shows that for our conditions the velocity is approximately proportional 

to the applied field. Tho E obtained by Karge, Hooks, and Oleson was 

^ max 

6 volts por on por mm Hg. This figure is well v/ithin the range of fields 

vdiorein exporimontal data reveals that velocity is proportional to field 

(fig. 3). No attempt has been made to justify or refute the m.f.p. 

( 4 ) 

difficulty. Von Engel' ‘ further states that the equation holds in the 
case of ions which movo and diffuse in moderately strong fioldsj in 
strong fields the velocity is proportional to tho square root of the field, 
lie fvirthor states that it can be assumed that the equation holds for weak 
fields. That our conditions fall within tho purview of Von Engel's state- 
ments can be seen from figure 4. wherein the velocity of tho ions is shovm 
to be proportional to the field. 

Graham and Rvihlig^^^^ obtained a value of 2 X 10® cm^/sec for the 
diffusion coefficient of the electron (D“) by using Droyvesteyn distri- 
bution in conjunction vfith the Bolt 2 mann transport equation. This valuo 
v/as calculated at the same conditions specified for our case. As can be 
seen from page 18 we calculated a value of 6.46 X 10^ cm^/sec. No 
attempt vj^as mado to resolve this difference. However, as vrill be shown 
later, this difference may bo rather easily invosi;igated by use of the 
program we propose by merely varying the coefficients of certain terms. 
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CliAPTI^iR IV 



PR0GRAI.!2'nNG 

In order to facilitate an understanding of the computer program we 
will set dov/n a brief mathematical description of what is talzing place 
within the computer* For convonience let if z and V = w-»-equations (4) 
nay then be written: 

i. = D-/^ - -[ h ■/ 4ireK7^(z - w)J z 

o X 

7^ / 2.345 X 10^^ 

= D"-^ -/K"(Eq ~ It ~ 4TreK“(z - w)J vr 

2.345 X 10^^ 

and from Poisson's 3qn, 

iii* - 41T e(z - vf) 

^x * " 

Furthermore, v/o arbitrarily establish as an initial condition: 

- 0, at X ts 0 

Let us no7)T divide the part of our positive column wherein our equations 
will apply into n equal parts* 

Let : 

Xq denote the origin of our column ( - O) 

Xj^ be the first point at a distance Ax from x^ 

then, 

Xj_ - iAx, where i n 0# 1# 

and, 

X ts Xj^ - Xj^_]L = ^ being the length of the column 

Let : 

Zi S z(Xi»t) 

Wi = w(xi,t) 

“ 4?Te(zi - (from iii.) 
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Romombering that it vms proposed to replace the partial derivatives with 
total derivatives (page 21 )* wo now observe the technique of doing this 
by the finite approximation formulae: 



^ u. / - u. 

iu^ _ \ • irl 1 



d w / \ 



AX 



\ 2 u. . — 2u . ^ u . , 

d u / \ . i /1 1 ~ "^i-l 

= x^) i i- 



^ x^ 



(ax)' 



Utilizing this expediency on our first three equations in this chapter, 
we have: 



i.a) 



ii*a) 



4j/eK'^(zj^ - VTj^^ / y,,w^/2.345 X 10^^ 

dw. -p|- 

— i _ - - (vr / - 2 w. T^w. 1 ) •/ (E •/ e^ )iv / - w. ) 

dt (ax) 2 irl ^ ^ A x' o ^ ^i ''i-/l i' 



-[y, - 4rreK"(z. - w.)jw. ■/ y^z. 7 ^ 2.345 X 



10 



14 



iii.a) e„ = e„ 7 ^ 4rre (^x) (z^ - w^ ) 

It is proposed to examine our solution over ten inoremonts; therefore, 
we restrict our problem as follows: 

i ** 1, 2, 3, •••• f 10 

Wii= zii Z 0 

Rearranging equations i.a), ii.a), and iii.a) above, we have: 
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- [^ 3 / 4jreK7^(zi - vr^)]J z^ 



2F7^ 






ii«b) 



dwj 

dT 



[f^ h-l >f "1 -^2.345 X 10 “ 



D- 

(Ax)‘ 






^ °xi> V [l - - ^^)JJ 

i-1 •/ 2.345 X 10^^ 



w. 



•/ 






iii.b) 



=»X 4 - 



(iix)2 

©x^_^ -f 4jTe(Ax)(zi - Wi) 



Let us now introduce some new constants by oombinins some of those which 

frequently occur; specifically lot, 

Vrf 

-g s .387 (arbitrarily setting x a 3cm) 



Ai- 



(ax) 



r?l= 1 . 



.2 - _ , ,.S5 X 10- 

A3 r 4T^eK'/ = 2.44 X 10"^ 
D“ 



Bl = 



(ax) 



^ 3 .720 X 10^ 



B2 : = .6 X 107 



B3 s 4ireK" s .109 

Inserting these new constants into equations i.b), ii.b), and iii.b) 
yields: 

- -^ 2(^0 ■/®Xi)Jzi/l 

7 ^®Xi) - 2A^ - A3jn; 7^ (Zi - w.)JJ z^ 

T^lH-l V>i 7^ 2.345 X 10^^ 
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uw . p 

ii.o) ~ a - ^ 2(^0 

“{^ 2(20 4- Qx^) -/2 Bi -/ Bgj^n- - ( 2 ^ - Wi)]jwi 
/ Bivri_i VBgnoZi -/ 2.345 X 10^^ 

iii.o) , •/ 6.023 X 10 "^(ax)(z. - w.) 

It should be noted that the value of n~ is arbitrarily takon to be 7 x IO^^q 
IT hile this may or may not be exact, it is the cor root order of magnitude# 

Further, in the ambipolar region s • In order to avoid overoomplioating 

the manipulation of decimals a new variable is introduced: 

Zi : lolSr; 

‘ Wj^ 5 

Changing variables and inserting the values of the constants vjhere applicable; 

^ \ - r 397' . /3 5 X (. + |. /35 0? ^ 6:t') 7?f-.2H X ?/ 

dt ^ ^ ^ 

f.387 t.l70S t.23i5 

^ t 

^ z[,7Z0 KIO'*-. 4X/F6 02 + e?;y)Jw.^ +{U/dV.02/-ex;) + /.V<f//0V,/09 m'%niO’t (i'- -W/jjJ w! 

+ .7-?OX/oV' +. 7i3 £'. t. 2345 

hi t 

fix-- ^ ,180? XI0^(2]^W.) 
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Lot us define s *02 7 ^ and S 2 . = then tho preceding 

equations can be v/ritten as: 

^ 1 3 5 X 5 , ; j ^.135 X Kfsri - . 77V - . V X /O'^/; 7 X/0^+ 5^ ;]J Z; 

f .I?d8 ^10* vJ^. -h , ZJ^S 

J!''' -- [, yZOK/o'*. .klMb* 5< , (.UKps,i f I0<f X (b‘^[, 7 X/^’^-5? i7} 

+ . 7-?0 X /£>" w:., t X /^'“h; t, ^5V5 






These are the final forms of the equations which the ocanputer %vill be 
asked to solve. A general idea of the method of solution oan be obtained 
by considering these equations in conjunction with the flov: diagram which 
appears on the foliating page. The details of execution of tho program 
are found follov/ing the flov^ diagram. 
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6;. . -- 0 


1 


Zi — W\ - S^i 


. 1807 X !0*Sz\ 


^.iu7 /lo^Szi ~ 


■ oi*e^c- S,; 


./35 X/0* Sn 


. . 339- ./as X/o^Sii 


a,-- 6389*./35 X/0^5n)z',v, 


.7XW'"' 


-.Z44 X/0'^(.7//0*V Sz.') 


-.77'/ - [:24-7X/0'^(.7X/0'V Szi)] 


J35 niO^S,i-(.l78i-C244 xiO'^r Szi)]} 


a2-C/35//0^S,i- x/O'USzi)!} ?')•*/ 


3j = .3?7 X 


3y - . 1708 y /o^ w,' 


a, t 


dj t ( a, ^ S 2 ) 


(5j -r 3^)7 


f [ay ^£aj + (a, + a2)]j 


. 6 X /o’ S/i' 


. 720 x/o"'- .tx/o’ 


6, - C 720 iC)’ 5/ { ) w}, , 


’7x/0'‘^- Szi 


J09 )(io'U.iy/0'‘'- Szl) 


/. vy x/oV C /07 X /o'®r,7 x/o*’- S 2 , )7 


,iX/0’’ Slit [i.n )i.ioUiioUio'U-7 /iO"* -Sii)] 


X [t X/o’S/i' tf/.VVX/0''+ [lOf X/fl'Y- 7X/i>*’'- Sz/)]}j 


- .720 X /O^ / w'/,, 


6y= .743 x/o'" 


b,*’ hz 


bj + (6,<- b,) 


by i-Cbz + ib, *bi)] 


,Z3¥£ +[bj^(b,i-bz)]} 


iH-^7 


irest my-/ > i 


,, /70 


A/uMe/*/ca? f fj t^o/? rouCi/)e (liongc* mctkoJ) 
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The following is the detailed program of the ccanputor solution to our 
problem. 



Loo 


Opn 


Opn 

Code 


M 

1 


M 

2 


M 

3 


Remarks 


0010 


md 


26 


2100 


2100 


0205 




0011 


ad 


35 


2100 


0201 


0200 


1 X 2"12 ^ 1 X 2-24 / i X 


0012 


ad 


35 


0200 


0200 


0207 




0013 


ad 


35 


0210 


0207 


0014 




0014 


su 


36 


[4 


1 

■wi 


020^ 


T t 

zi - MTi s in 0202 


0015 


tm 


34 


0202 


2100 


0020 


no -*0016, yes -*•0020' 


0016 


md 


26 


2100 


2100 


2000 


n _(2000 

(2001 


0017 


tm 


34 


3000 


2100 


0022 




0020 


sf 


31 


0240 


0202 


20001 


.1807 X 10® X Sji in [goS 


0021 


mr 


25 


0241 


2001 


2001 J 


0022 


ad 


35 


2100 


0205 

0206 


2002^ 

2003 / 


Siooi 


0023 


ad 


35 


2100 


0024 


tm 


34 


3000 


2100 


0400 


gives in [2000 


0025 


ad 


35 


2100 


0242 


2002A 


[lo03 


0026 


ad 


35 


2100 


0243 


2005 J 


0027 


tm 


34 


3000 


2100 


0400 


.02 -/e^^ r in [|°°J 


0030 


ad 


35 


2100 


2000 


0203-\ 


in 0203 and 0204 


0031 


ad 


35 


2100 


2001 


0204 j 


0032 


tm 


34 


2001 


2100 


0035 


no -t 0033, yes 0035 


0033 


md 


26 


2100 


2100 


2000 


0 2000, 2001 

% 


0034 


tm 


34 


3000 


2100 


0037 




0035 


ad 


35 


2000 


0244 


2000 ^ 












0245 




-.135 X 10° 5 ^^ 


0036 


mr 


25 


2001 


2001 J 





-24 
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0037 


ad 


35 


2100 


0246 


2002 




0040 


ad 


35 


2100 


0247 


2003 


•389 in 2002 and 2003 


0041 


tm 


34 


3000 


2100 


0400 


.389 - .135 X 10® sii in [goS 


0042 


tm 


34 


2001 


2100 


0045 


no -»0043, yos^-0045 


0043 


md 


26 


2100 


2100 


2000 




0044 


tm 


34 


3000 


2100 


0054 




0045 


ex 


32 


0014 


0211 


0212 




0046 


ad 


36 


0212 


0213 


2004 




0047 


ex 


32 


2004 


0211 


0050 




0050 


tm 


34 


[kA] 


2100 


0053 


no -*0061, yes-*0053 


0051 


md 


26 


2100 


2100 


2000 




0052 


tm 


34 


3000 


2100 


0054 




0053 


tm 


34 


3000 


2100 


0500 


a, in (2000 
®1 ^ (2001 


0054 


ad 


35 


2100 


2000 


0214 


an in 0214 
1 (0215 


0055 


ad 


36 


2100 


2001 


0215 




0056 


ad 


35 


2100 


2100 


2000'] 


^ in 2000 and 2001 


0057 


ad 


35 


2100 


0202 


200lJ 


0060 


ad 


36 


2100 


0250 


20021 


1- .7 X 10"^ in 2002 and 2003 


0061 


ad 


35 


2100 


0251 


2003 J 




0062 


tm 


34 


3000 


2100 


0400 


.7 X 10"“^ / Bji in 2000 and 2001 


0063 


tm 


34 


2000 


2100 


0066 


no -*0064, yes -* 0066 


0064 


md 


26 


2100 


2100 


2000 




0065 


tm 


34 


3000 


2100 


0072 




0066 


ad 


35 


0252 


2000 


2000' 


1 .244 X 10^2 (^7 X 10"“^ / 6^^) 


0067 


mr 


25 


0253 


2001 


2001, 


0070 


ad 


35 


2100 


0254 


2002^ 


.774 in 2002 and 2003 


0071 


ad 


35 


2100 


0255 


2003 J 


1 



38 




I 

I 



I 

I 

I 



I 



I 






I 



0072 


tm 


34 


3000 


2100 


0400 


-.774 -.244 X 10^2 (^7 ^ 10"^ / s^i) in 


0073 


ex 


32 


0212 


0211 


0074 


2000 and 2001 


0074 


tm 


34 


f'l] 


2100 


0077 


no -^0075, yes — *• 0077 


0076 


md 


26 


2100 


2100 


2000 




0076 


tm 


34 


3000 


2100 


0101 




0077 


tm 


34 


3000 


2100 


0510 




0100 

0101 


ad 


35 


2100 


2000 


0216^ 


, (0216 
^2 “ (0217 


0102 


ad 


35 


2100 


2001 


0217/ 


0103 


ad 


35 


2100 


0257 


2001 


.387 in 2001 


0104 


ad 


35 


2100 


0212 


2004 




0105 


su 


36 


2004 


0213 


2004 




0106 


ox 


32 


2004 


0211 


0107 




0107 


tm 


34 1 


fii-i] 


1 2100 


0112 


no -^OllO^ yes -*'0112 


0110 


md^ 


26 


2100 


2100 


2000 




0111 


tm 


34 


3000 


2100 


0114 




0112 


tm 


34 


3000 


2100 


0516 




0113 

0114 


ad 


35 


2100 


2100 


0220 


a3 in 0220 and 0221 


0115 


ad 


35 


2100 


2001 


0221 


a3 in 0220 and 0221 


0116 


ad 


35 


0212 


0226 


2004 




0117 


ex 


32 


2004 


0211 


0120 




0120 


tm 


34 ‘ 


[wl; 


j 2100 


0123 


no -► 0121, yes ^0123 


0121 


md 


26 


2100 


2100 


2000 




0122 


tm 


34 


3000 


2100 


0125 




0123 


tm 


34 


3000 


2100 


0524 





39 



t 



\ 




0124 



0125 


ad 


35 


2100 


2000 


0222" 


I B4 in 0222, 0223 


0126 


ad 


35 


2100 


2001 


0223. 


0127 


ad 


35 


2100 


0220 


2002'j 


f ag in 2002, 2003 


0130 


ad 


35 


2100 


0221 


2003 i 




0131 


tm 


34 


3000 


2100 


0400 


B3 / 04 in 2000, 2001 


0132 


ad 


35 


2100 


0216 


2002 


02 in 2002, 2003 


0133 


ad 


35 


2100 


0217 


2003 


0134 


tm 


34 


3000 


2100 


0400 


^■2 7^ ^3 7^ ®^4 2000, 2001 


0135 


ad 


35 


2100 


0214 


2002" 


i in 2002, 2003 


0136 


ad 


35 


2100 


0215 


2003 J 


0137 


tm 


34 


3000 


2100 


0400 


^ 2 . ^ 2000^ 2001 


0140 


ad 


35 


2100 


0262 


2002' 


1 .2345 in 2002, 2003 


0141 


ad 


35 


2100 


0263 


2003 J 




0142 


tm 


34 


3000 


2100 


0400 


^ in 2000, 2001 
dt 


0143 


ad 


35 


2000 


0227 


2000 


^ scaled) in 2000, 2001 


0144 


tm 


34 


3000 


2100 


0145 




0145 


ta 


33 


0157 


2000 


0150 




0146 


pr 


21 


2000 


0147 


0002 




0147 


ht 


22 


f 


f 


f 




0150 


sm 


30 


2001 


2000 


2002 


in 2002 

^ Z4 J in m3 of 2004 


0151 


sm 


30 


0014 


0227 


2004 


0152 


ox 


32 


2004 


0156 


0154 




0153 


ad 


35 


0154 


0277 


0154 


in m3 of 0122 


0154 


ad 


35 


2002 


2100 


f 


dz^* 

transferred to integration sub-routine 


0155 


tm 


34 


3000 


2100 


0160 





40 




I 



I 



I 



vl 



0156 






f 


f 


7777 


0157 






f 


f 


0001 


0160 


tm 


34 


0204 


2100 


0163 


0161 


md 


26 


2100 


2100 


2000 


0162 


tm 


34 


3000 


2100 


0170 


0163 


ad 


36 


0203 


0264 


2000 


0164 


mr 


25 


0204 


0265 


2001 


0165 


su 


36 


2100 


2001 


2001. 


0166 


ad 


36 


2100 


0266 


2002 


0167 


ad 


36 


2100 


0267 


2003 


0170 


tm 


34 


3000 


2100 


0400 


0171 


tm 


34 


2001 


2100 


0174 


0172 


md 


26 


2100 


2100 


2000 


0173 


tm 


34 


3000 


2100 


0311 


0174 


ad 


35 


0120 


0213 


2004 


0175 


ex 


32 


2004 


0211 


0304 


0176 


tm 


34 


3000 


2100 


0304 


0177 












0200 












0201 






0001 


0001 


f 


0202 












0203 












0204 












0205 












0206 












0207 













no -►0161, yes-rOISS 



.6 X lO' in 2000, 2001 

.720 X 10^ in 2002, 2003 

.720 X 10“^ -.6 X lO"^ sii in 2000, 
no -»0172, yes -►0174 



i X 2’12 yl i X 2“24 



temp storage 



} 



'li 






2 X (i X 2-12 y: i X 2“^^) 



2001 



41 



0210 

0211 

0212 

0213 

0214 

0216 

0216 

0217 

0220 

0221 

0222 

0223 

0224 

0226 

0226 

0227 

0230 

0231 

0232 

0233 

0234 

0236 

0236 

0237 

0240 

0241 



36 



26 



1305 


1306 


0202 


dummy for 


0014 


7777 


f 


f 


extracted 


for 0034 


f 


2001 


2001 


diJimniy for 


0036 


0002 


f 


± 







0001 
02 f 



J a-j^ temporary storage 

^ a£ temporary storage 

^ ag temporary storage 

I’ temporary storage 

L temporary storage 

1 ^ 



f f 
f. 0030 



f f 0033 ' 

1047 3476 0000 > 



storage 
bg storage 
bg storage 
b^ storage 
.1807 X 10^ 



42 



0242 

0243 

0244 

0245 

0246 

0247 

0250 

0261 

0252 

0253 

0254 

0255 

0256 

0257 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0270 

0271 

0272 

0273 



02 



02 



02 



02 



f 


f 


f 


0121 


7300 


0000 


f 


f 


0017 


3227 


3777 


7720 


f 


f 


f 


3071 


2601 


0142 


f 


f 


f 


0000 


2226 


3173 


f 


f 


0047 


3441 


3567 


6333 


f 


f 


f 


6142 


2335 


1362 


f 


f 


f 


3061 


1156 


4571 


f 


f 


0033 


1011 


1730 


0000 


f 


X 


f 


1700 


4061 


1156 


f 


f 


0030 


2670 


6600 


0000 


0000 


0000 


0017 


1604 


0000 


0000 


f 


f 


0060 


3062 


1212 


2176 


f 


f 


0017 


3407 


7777 


4512 



-.135 X 10^ 
.389 

.7 X 10”^ 
-2.44 X 10^2 
-.774 
.387 

.1708 X 10® 
.2345 
.6 X lo"^ 
.720 X 10^ 
.109 X 10^5 
1.44 X 10^ 



43 



0274 

0275 

0276 

0277 

0300 

0501 

0302 

0303 

0304 

0305 

0306 

0307 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0520 

0321 

0322 

0323 

0324 

0325 

0326 



f 

7066 


f 

2107 


0041') 
7440 j 


.763 X 


10^° 


0013 


0013 


f 


m -/ 1 


m /l f, m - 12 foctij) 


f 


f 


0027 


(2^7^!) 


X 2"36 



■bm 


34 


J 2100 


0307 


no -► 0305, yBS -^0307 


md 


26 


2100 2100 


2000 




'tan 


34 


3000 2100 


0311 




tm 


34 


3000 2100 


0540 





ad 


36 


2100 


2000 


0230') 












C b-, in 0230, 0231 




ad 


35 


2100 


2001 


0231 j 




su 


36 


2100 


0202 


2001 Sjjj^ in 2001 




ad 


35 


2100 


2100 


2000 




ad 


35 


2100 


0250 


2002") 












\ .7 X 10“^ in 2002, 2003 




ad 


35 


2100 


0251 


2003 ) 




tm 


34 


3000 


2100 


0400 * 7 X 10*"^ - s^i in 2000, 


2001 


tm 


34 


2001 


2100 


0323 no -*'321, yes ■♦323 




md 


26 


2100 


2100 


2000 




tm 


34 


3000 


2100 


0330 




ad 


35 


2000 


0270 














1C A 

1 .109 X 10-*-° X (.7 X 10"^ 


- Sj,i) 


rar 


25 


2001 


0271 


2001 J 




ad 


35 


2100 


0272 


2002") 












L lo44 X 10^ in 2002, 2003 




ad 


35 


2100 


0273 


2003 / 





44 



0327 


tm 


34 


3000 


2100 


0400 


1.44 X 10^7^ .109 X 10^5 X (.7 X 10"^-Sjji) 














in 2000, 20001 


0330 


tm 


34 


0204 


2100 


0333 


no-*-0331, yos-^-0333 


0331 


md 


26 


2100 


2100 


2002 




0332 


tm 


34 


3000 


2100 


0336 




0333 


ad 


35 


0203 


0264 


2002 




0334 


mr 


25 


0204 


0265 


2003 




0335 


tm 


34 


3000 


2100 


0400 




0336 


ex 


32 


0120 


0211 


0337 


1 


0337 


tm 


34 


["!] 


2100 


0342 


no 0340, yo s -»• 0342 


0340 


md 


26 


2100 


2100 


2000^ 




0341 


tm 


34 


3000 


2100 


0344 




0342 


tm 


34 


3000 


2100 


0550 


y 1>2 in 0252, 0233 


0343 














0344 


ad 


35 


2100 


2000 


0232 




0345 


ad 


36 


2100 


2001 


0233 j 




0346 


ad 


35 


2100 


0266 


2000 


i .720 X 10“^ in 2000, 2001 


0347 


ad 


35 


2100 


0267 


2001 J 




0350 


su 


36 


0120 


0213 


2004 




0551 


ex 


32 


2004 


0211 


0352 




0352 


tm 


34 


[wl-f 


1 2100 


0355 


no 0363, yos-^0355 


0353 


md 


26 


2100 


2100 


2000 




0354 


tm 


34 


3000 


2100 


0357 




0355 


tm 


34 


3000 


2100 


0560 




0356 














0357 


ad 


35 


2100 


2000 


0234' 


1 b3 in 0234, 0235 


0360 


ad 


35 


2100 


2001 


0235. 



45 



0361 


ad 


35 


2100 


0274 


2000 "i 


y .763 X 10^° in 2000, 2001 


0362 


ad 


36 


2100 


0275 


2001 - 


) 


0363 


ox 


32 


0074 


0211 


0364 




0364 


tm 


34 


[4] 


[ 2100 


0367 


no -*0365, yes -»0367 


0365 


md 


26 


2100 


2100 


2000 




0366 


tm 


34 


3000 


2100 


0371 




0367 


tm 


34 


3000 


2100 


0570 




0370 














0371 


ad 


35 


2100 


2000 


0236^1 


. 64 in 0236, 0237 


0372 


ad 


35 


2100 


2001 


0237 j 


0373 


ad 


35 


2100 


0234 


2002^ 


\ 63 in 2002, 2003 


0374 


ad 


36 


2100 


0235 


2003 J 


0375 


tan 


34 


3000 


2100 


0400 


1>3 7^ ^4 in 2000, 2001 


0376 


tm 


34 


3000 


2100 


0440 




0377 










« 




0400 


si 


27 


3000 


0422 


2004 




0401 


ad 


35 


0423 


2004 


0410 




0402 


tm 


34 


3000 


2100 


0426 




0403 


su 


36 


2002 


2000 


2004 




0404 


sm 


30 


2003 


2004 


2003 




0406 


ad 


35 


2003 


2001 


2001 




0406 


to 


37 


2001 


3000 


0416 




0407 


sf 


31 


2000 


2001 


2000 




0410 












exit 


0411 


su 


36 


2000 


2002 


2004 




0412 


sm 


30 


2001 


2004 


2001 





46 



0413 


ad 


36 


2002 


2100 


2000 






0414 


tm 


0 4 


3000 


2100 


0405 






0415 


si 


27 


2001 


0424 


2001 






0416 


sm 


30 


2001 


0424 


2001 






0417 


si 


27 


2001 


0425 


2001 






0420 


ad 


35 


0425 


2000 


2000 






0421 


tm 


34 


3000 


2100 


0407 






0422 




02 


f 


f 


0030 






0423 


tm 


34 


3000 


2100 


f 






0424 




02 - 


f 


f 


0001 






0425 






f 


f 


0001 






0426 


tm 


34 


2001 


2100 


0431 






0427 


tm 


34 


2003 


2100 


0433 






0430 


tm 


34 


3000 


2100 


0410 






0431 


sf 


31 


2000 


2001 


2000 






0432 


tm 


34 


3000 


2100 


0427 






0433 


sf 


31 


2002 


2003 


2002 






0434 


ta 


33 


2002 


2000 


0411 






0435 


tm 


34 


3000 


2100 


0403 






0436 
















0437 
















0440 


ad 


35 


2100 


0232 


20021 


b 2 in 2002, 


2003 


0441 


ad 


35 


2100 


0233 


2003 j 






0442 


tm 


34 


3000 


2100 


0400 


^2 7^ ^3 Ad 


in 2000, 2001 


0443 


ad 


35 


2100 


0230 


2002*) 


bi in 2002, 


2003 


0444 


ad 


35 


2100 


0231 


2005 J 







47 



0445 


tm 


34 


3000 


2100 


0400 


1>1 /■ t)2 •/ ^3 7^ ^4 2000, 2001 


0446 


ad 


35 


2100 


0262 


2002'^ 


,2345 in 2002, 2003 


0447 


ad 


35 


2100 


0263 


2003 J 




0450 


tm 


34 


3000 


2100 


0400 


t 

dw4 














in 2000, 2001 


0451 


ad 


35 


2000 


0227 


2000 


dw^ 
















0452 


ta 


33 


0157 


2000 


0455 




0453 


pr 


21 


2000 


0454 


0002 




0454 


ht 


22 


f 


f 


f 




0455 


sm 


30 


2001 


2000 


2002 




0456 


sm 


30 


0337 


0227 


2004 




0457 


ex 


32 


2004 


0156 


0461 




0460 


ad 


35- 


0461 


0277 


0461 




0461 


ad 


35 


2002 


2100 


f 


dVT* 

transferred to into, ’ration sub-routine 


0462 


ad 


35 


0201 


0200 


0200 


(i/-l)X 2“2'i^iX2'^V 














iX2-24 


0463 


tm 


34 


0276 


0200 


0012 


test m 7 ^ 1 > i 


0464 


tm 


34 


3000 


2100 


1100 


U.T. to integration sub-routine 


0465 


tm 


34 


3000 


2100 


0010 




0500 


sm 


30 


0050 


0506 


2005 


[_2iy'lJ in mg of 2005 


0501 


ex 


32 


2005 


0507 


0502 




0502 


sf 


31 


2100 


•H 


2002 




0503 


ad 


35 


2000 


2002 


2000 ^ 


a]_ in 2000, 2001 


0504 


mr 


25 


2001 


2003 


ZOOlJ 


0505 


tm 


34 


3000 


2100 


0054 




0506 




02 


f 


f , 


0014 




0507 






f 


7777 


f 




0510 


sm 


30 


0212 


0506 


2005 





48 



1 



I 



I 

I 



... 



0511 


OX 


. 32 


2005 


0507 


0512 








0512 


sf 


31 


2100 [ 




2002 








0513 


ad 


35 


2000 


2002 


2000 \ 


®2 


in 


2000, 


0514 


inr 


25 


2001 


2003 


2001 j 






0515 


tm 


34 


3000 


2100 


0101 








0516 


sm 


30 


0107 


0506 


2005 








0517 


ex 


32 


2005 


0507 


0520 








0520 


sf 


31 


2100 j 


■ 1 ~ 


1 2002 








0521 


ad 


35 


2000 


2002 


2000 'j 


»3 


in 


2000, 


0522 


znr 


25 


2001 


2003 


2001 / 






0523 


tm 


34 


3000 


2100 


0114 








0524 


sm 


30 


0120 


0506 


2005 








0525 


ex 


32 


2005 


0507 


0526 








0526 


sf 


31 


2100 




2002 








0527 


ad 


35 


2000 


2002 


2000 'j 


&4 


in 


2000. 


0530 


mr 


25 


2001 


2003 


2001 J 






0531 


tm 


34 


3000 


2100 


0125 








(cells 0532 through 0537 are filled with 


L zeros) 


0540 


sm 


30 


0304 


0506 


2005 








0541 


ox 


32 


2005 


0507 


0542 








0542 


sf 


31 


2100 1 


iw[^ 2002 








0543 


ad 


35 


2000 


2002 


2000'\ 


^1 


in 


2000 


0544 


mr 


25 


2001 


2003 


2001} 






0545 


tm 


34 


3000 


2100 


0311 









0546 

0547 



2001 



2001 



2001 



2001 



49 




1 

i 

I 

I 

I 



I 



0550 


sm 


30 


0337 


0506 


0551 


ex 


32 


2005 


0507 


0552 


sf 


31 


2100 


r™i] 


0553 


ad 


35 


2000 


2002 


0554 


mr 


25 


2001 


2003 


0555 


tm 


34 


3000 


2100 


0556 










0557 










0560 


rrm 


30 


0352 


0506 


0661 


ex 


32 


2005 


0507 


0562 


sf 


31 


2100 


'^i-l] 


0563 


ad 


35 


2000 


2002 


0564 


mr 


25 


2001 


2003 


0565 


tm 


34 


3000 


2100 


0566 










0567 










0570 


sm 


30 


0364 


0506 


0571 


ex 


32 


2005 


0507 


0572 


sf 


31 


2100 


[h] 


0573 


ad 


35 


2000 


2002 


0574 


mr 


25 


2001 


2003 


0575 


tm 


34 


3000 


2100 



0576 

0577 



2005 

0552 

2002 

2000 

i bp in 2000, 2001 

2001 j 

0344 

2005 

0562 

2002 
2000 ") 

I b 3 in 2000, 2001 
2001 J 

0357 

2005 

0572 

2002 

2000 

2001 

0371 



50 



1000 



1001 

1002 

1003 

1004 

1005 

1006 
1007 
1010 
1011 
1212 

1013 

1014 

1015 

1016 
1017 
1020 
1021 



sl 


27 


3000 


1016 


2001 


ad 


35 


1014 


2001 


1013 


od 


35 


2000 


2100 


1022 


bl 


05 


3000 


3000 


1014 


ex 


32 


2002 


1000 


2000 


mr 


25 


2002 


2007 


2002 


ex 


32 


2002 


2001 


2000 


ex 


32 


2100 


2001 


2002 


sl 


27 


2001 


2005 


2001 


tm 


34 


2001 


2100 


1005 


pr 


21 


2000 


1020 


0001 


tm 


34 


3000 


2100 


f 




02 


f 


f 


0004 




02 


f 


f 


0030 






5000 


f 


f 




01 


f 


f 


f 






7400 


f 


£ 



(cells 1022 through 1047 are filled with zeros) 



1050 


ad 


35 


2100 


2100 


1057 


clear i tally 


1051 


ad 


35 


1060 


1057 


1052 




1052 


^ad 


35 


1305 /i 2100 


2000j 




1053 


tm 


34 


3000 


2100 


1000 


U.T. to convert and print sub-routine 


1054 


ad 


35 


1062 


1057 


1057 


i ■/ 1 -♦ i 


1055 


tm 


34 


1061 


1057 


1051 




1056 


tm 


34 


3000 


2100 


1160 





51 




I 

I 



I 

I 

I 



1057 



i X 2-12 



1060 


ad 


35 


1305 


2100 


2000 




1061 






0030 


. f 


f 


m 7 ^ 1 ff 


1062 






0001 


f 


f 




(Tho follOTvlng 


is the 


sub-routine for 


the integration of n 1 st order 


equations 


by the Runge 


-Kutta-Gill method) 


1100 


Si 


27 


3000 


1242 


20001 


Prepare exit to main program 


1101 


ad 


35 


1246 


2000 


1160 J 




1102 


tm 


34 


1163 


1303 


1200 


Test if first time through 


1103 


tm 


34 


3000 


2100 


1107 


U. T. 


1104 


ad 


35 


2100 


2100 


1247 


0 -»'k (clear print tally) 


1105 


ad 


35 


2100 


2100 


1250 


0 -#► j (clear iteration tally) 


1106 


ad 


35 


2100 


1177 


1251 


1 X 2“^^ (reset var. no* tally) 


1107 


ad 


35 


1250 


1251 


1252 


i X 2-12 ^ j X 2-24 tally 


1110 


ad 


35 


1234 


1262 


1111 


[ 9 J ] in 11 and in of 1111 


1111 


l^sm 


30 


4 


bJ 


200l] 




1112 


ad 


35 


1232 


1251 


1113 


ik\] in Ml of 1113 


1113 j 


[su 


36 




2001 


2 OO 4 J 




1114 


ad 


35 


1267 


1250 


1115 


[A^ ] in II 2 of 1116 


1115 ^ 


l^mr 


26 


2004 


Ad 


2005] 




1116 


su 


36 


2004 


2005 


2006 




1117 


to 


37 


2006 


2100 


1226 


test (truncation tost) 


1120 


mr 


25 


2006 


1301 


2007 




1121 


sm 


30 


1251 


1242 


1270 


i X 2"^® from i X 2“^^ 


1122 


ad 


35 


1251 


1270 


1271 


i X 2“^2 ^ i X 2“^® 


1123 


ad 


35 


1272 


1271 


1124 


J in and Kg of 1124 



52 




I 



1124 


r ad 


35 




2007 


xJ 1 


xj X AXJ - XJ “/ 1 X j 




L 




1 




} ' 


i i - i i 


1125 


sm 


30 


1251 


1235 


2003 


i X 2-24 from i X 2-12 


1126 


ad 


35 


1251 


2003 


1273 


i X yf i X 2-24 


1127 


ad 


35 


1237 


1273 


1130 


[q'!] ,in M]_ and [^■‘J in Mg of 1130 


1130 


/su 


36 


■q^ 


• . 


20071 






L 




1 


1 


J 




1131 


ad 


35 


1241 


1252 


1132 


[X- j In Ml and £c'J in Mg of 1132 


1132 




25 


1 J 
1 


cj 


2000 j 


• 4 

X. X d 


1133 


sm 


30 


2007 


1274 


2007 


/cjf-xj) X 8-1 


1134 


sm 


30 


2000 


1274 


2000 


xl xd X 8 “^ 


1135 


ad 


35 


2000 


2007 


2001 


[ijKc'-O iU X8-1 


1136 


mr 


25 


2006 


1275 


2006 


3 r/ X 8 "^ 


1137 


ad 


35 


2001 


2006 


2002 




1140 


sm 


30 


1111 


1242 


2003 


in M 3 of 2003 


1141 


ad 


35 


1276 


2003 


1142 


in M 3 of 1142 


1142 




30 


2002 


1277 


?i]- 




1143 


ad 


35 


1177 


1251 


1251 


(i/l) X i X rl2 


1144 


tm 


34 


1300 


1251 


1107 


test (n 7 ^ 1 ) >1 


1145 


ad 


35 


2100 


1177 


1251 


1 X 12"^2 i X 2"^2 


1146 


ts 


17 


2010 


3000 


1164 


print j X 2”24 if ts 2010 up 


1147 


ad 


35 


1250 


1175 


1250 


(i/l) X X 2’24 


1150 


tm 


34 


1176 


1250 


1160 


test 4 > j 


1151 


ad 


35 


2100 


2100 


1250 


0 ^ j X 2"24 


1152 


ts 


17 


2020 


3000 


1166 


accumulate X if 2020 up 


1153 


ad 


35 


1163 


1247 


1247 


k-/ 1 k 


1154 


tm 


34 


1302 


1247 


1160 


test K ^k 


1155 


ad 


35 


2100 


2100 


1247 


0 -»-k 


1156 


ts 


17 


2040 


3000 


1170 


print Xf octal if 2040 up 
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I 



1 



I 



1157 



1160 

1161 

1162 

1163 

1164 

1165 

1166 
1167 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 
1200 
1201 
1202 

1203 

1204 

1205 

1206 
1207 
1210 



tm 


34 


3000 


2100 


1050 


U, T, to nev/ print routine 












U. T. to SR 






f 


f 


0001 


1 X 2“^® 


pr 


21 


1250 


2100 


0001 


print 3 X 2“^'^ 


tm 


34 


3000 


2100 


1147 


U. T. 


ad 


35 


1301 


1304 


1304 


X^/ (m/Df-. Xq-/ t 


tm 


34 


3000 


2100 


1153 


U. T. 


ts 


17 


2020 


3000 


1173 




pr 


21 


1305 


2100 




print Xj^, Xg^ .... ^X^ 


tm 


34 


3000 


2100 


1160 




pr 


21 


1304 


2100 


n/l] 


print Xq - t, X]_. X£, 


tm 


34 


3000 


2100 


1160 
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1 X 2”24 
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4 X 2“24 






0001 


f 


f 
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ad 


35 


1230 


1300 


2000 


fXo] in Hi of 2000 


ad 


35 


1231 


2000 


1232 


^'^1 ol* nod. dvnnray 1 


ad 


35 


1300 


2000 


2002 


in Ml of 2002 


ad 


35 


1233 


2002 


1^4 


J in Mj of nod. dummy 2 


sia 


30 


2000 


1235 


2003 


in Mg of 2003 


ad 


35 


2002 


2003 


2004 


[^^£>1 ^-2 2004 


ad 


35 


1236 


2004 


1237 


^'"1 nod.dicnmy 3 


ad 


35 


1240 


2000 


1241 


ill % oi* nod, dianny 4 


sm 


30 


2002 


1242 


2002 


J in M 3 of 2002 
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I 



1211 



sm 30 1300 1242 2000 n V- 1 X 2“^® in 2000 



1212 


ad 


35 


1243 


2000 


1173 


sets up print of X^, 


1213 


su 


36 


2000 


1163 


2001 


n -X 2 "^® in 2001 


1214 


ad 


35 


1244 


2001 


1171 


sets up print of , .... 


1215 


ad 


35 


2100 


2100 


1270 


0 i X 2 “^® 


1216 


ad 


35 


2002 


1270 


2005 


[<=i;J in I.I 3 of 2005 


1217 


ad 


35 


1245 


2005 


1220 


f«l . j in M 3 of 1220 


1220 


|ad 


35 


2100 


2100 




0 -qi 


1221 


ad 


35 


1163 


1270 


1270 


(i/-l) X 2"36 -)► 1 X 2"^® 


1222 


tm 


34 


1270 


2000 


1224 


tost i ^ n /l 


1223 


tm 


34 


3000 


2100 


1216 




1224 


ad 


35 


2100 


1163 


1303 


signal for first timo through 


1225 


tm 


34 


3000 


2100 


1104 


U. T, 


1226 


pr 


21 


2006 


1227 


0001 


print r^ 


1227 


ht 


22 


f 


f 


f. 


halt 


1230 






1304 


f 


f 




1231 


su 


36 


f 


2001 


2004 


dummy 1 


1232 ( 


su 


36 




2001 


2OO4J 


modified dummy 1 for 1112 


1233 


sm 


30 


f 


1253 


2001 


dxanmy 2 


1234 


^sm. 


30 


^lo 


1253 


2001 ] 


modified dummy 2 for 1111 


1235 




02 
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f 


0014 
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su 
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2007 


dummy 3 


1237 , 


[^su 


36 


qo 




2007 J 


modified dummy 3 for 1130 
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mr 
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2000 
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mr 


25 


• 


1263 


2OO0J 


modified dummy 4 


1242 




02 
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f 


0030 
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1243 

1244 

1245 

1246 
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1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 

1260 

1261 

1262 

1263 

1264 

1265 

1266 

1267 

1270 

1271 

1272 

1273 

* 

1274 

1275 



pr 


21 


1304 


2100 


f 


dxnnmy print for 1173 


pr 


21 


1305 


2100 


f 


dxjinniy print for 1171 


ad 


35 


2100 


2100 


f 


dxjinniy for 1220 


tm 


34 


3000 


2100 


f 


dummy exit oamaand 



k X 2"^® (print tally) 
j X 2“^^ (iteration tally) 
i X 2“12 (T^ar. no. tally) 
i X 2"^2 ^ j X 2"24 
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0001 
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b1 
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b2 
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0001 


b3 




4000 
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f 


1 — 1 

It 

o 

< 




5520 


2357 


4151 


a} = YiTz 


02 


5520 


2357 


4151 


a2 : 




6525 


2525 


2523 


= 5/6 




4000 


f 


f 


c° - 1/2 




5520 


2357 


4151 


It 

1 — { 

0 


02 


5620 


2367 


4151 


it 
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4000 


f 


f 


C3 = 1/2 


mr 26 


2004 


1257 


2005 


dianmj' for 1116 
i 3E 2”^S 

i X 2"^2 7^ i X 2"2‘ 


ad 35 


1304 


2007 


1304 


dvanmy for 1124 
i X 2"^^ / i X 


02 


f 


f 


0003 






3000 


f 


f 


3 X 8*^ 



56 



1276 sm 30 2002 1277 £ 

1277 f f 0003 

(Colls 1300 through 1332 were filled by tape. These cells contain the 
values of t, etc. These values are not included here because they 

depend upon the values of computed paraiaeters~see Chapter V») 

Seme notes on the program: 

All addresses appearing in brackets are not put in by punch cards but 
rather by other parts of the program. They are included for purposes of 
clarity in deciphering the program. 

The following notes apply to the integration sub-routine; 
ts 2010 — up for iteration number tally 
ts 2020 — up for generation of 
ts 2040 ~ up for print of x^ octal 
n ■ no. of dependent variables (equatiens) 

TT" = incremental change ef independent variable 
K = readout at every K point 
X|^“ independent variable 

ilexT print or readout routine to start at 1050. 

Automatic print of cell 2006 indicates interval of integration too 
large. 

Error of integration is proportional to the 5th povror of t. 
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ClIAPTiiR V 



CONPLUDING REUARKS AlID RECOlB.IBHDATIOilS 
FOR FURTHER WORK 

As noted previously, it vrould bo fatuous at this point to draw any 
definite oonolusions suient tho solution to oxur problem proposed in tho 
first chapter* On the basis of the ccanputer solutions we cannot say with 
absolute certainty that tho results are compatible with experimental 
observations, Vfe can say, howevor, on the basis of our sketchy results 
that tho theory and method are promising, Vfithout equivocation we can 
say that the computer program is eminontly successful— it does work and it 
does solve the equations proposed as a description of tho traveling 
density waves, A great deal of work is forecast in order to establish 
the identity of those equations vrith the phenomenon in the positive column. 
Great numbers of initial conditions should be inserted in the program and 
the o*oefficionts should be veiried in order to eventually establish the 
validity of the equations. It is believed that if this is done, the 
resulting solutions by the oomputor will correspond to the observations 
made under actual conditions in the laboratory, IlmTever, no precise 
conclusions can be drawn from a very few computations. Those computations 
cannot be made in a few hours. It is predicted that tho collection of 
enough data v/ill require on the order of several hundred hovirs of computer 
time. The important point to bo remembered is that a workable program is 
available and patience and perseverance will revj’ard the eventual researcher 
v/ith tho answer. It appears from initial computer runs that tho term 
zUq is dominated in the equations (4) by other terms v/hen X>*s n^. 

The initial computer runs suggest that "floating point arithmetic" might 
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I 



I 



bo necessary throughout the program duo to the drastic changes noted in 
derivatives and the time required for oomputation using the present 
program. Ho other trends vrere noted from the initial oomputor runs. 
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